Layered ferromagnet-superconductor structures: the tt state and proximity effects 
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We investigate clean mutilayered structures of the SFS and SFSFS type, (where the S layer is 
intrinsically superconducting and the F layer is ferromagnetic) through numerical solution of the 
self-consistent Bogoliubov-de Gennes equations for these systems. We obtain results for the pair 
amplitude, the local density of states, and the local magnetic moment. We find that as a function 
of the thickness dp of the magnetic layers separating adjacent superconductors, the ground state 
energy varies periodically between two stable states. The first state is an ordinary "0-state", in 
which the order parameter has a phase difference of zero between consecutive S layers, and the 
second is a "7r-state", where the sign alternates, corresponding to a phase difference of n between 
adjacent S layers. This behavior can be understood from simple arguments. The density of states 
and the local magnetic moment reflect also this periodicity. 

PACS numbers: 74.50+r, 74.25.Fy, 74.80.Fp 



I. INTRODUCTION 

The study of layered ferromagnet-superconductor 
(F/S) heterostructure has sustained the active interest of 
many researchers. This is due in great part to continuing 
and recent progress in the preparation and fabrication of 
multilayer systems, and to the potential use of such het- 
erostructures in various important applications. In par- 
ticular, structures consisting of alternating ferromagnet 
(F) and superconductor (S) layers may exhibit, in certain 
cases, a ground state in which the difference A</> between 
the order parameter phase of adjacent superconductor 
layers equals n. These are the so called u n junctions". 
These F/S hybrid structures offer advances in the field 
of nanoscale technology, including quantum computing, 1 
where the implementation of a quantum two-level sys- 
tem is based on superconducting loops of tt junctions. 
Furthermore, artificial composites involving a supercon- 
ductor sandwiched between two ferromagnets, the design 
of which follows from giant magnetoresistive (GMR) de- 
vices, show potential use as spin- valves^ and nonvolatile 
memory elements^ An essential principle behind many 
of these spin-based devices is the damped oscillatory na- 
ture of the Cooper pairs in the ferromagnet region, and 
the associated phase shift in the superconducting order 
parameter. 

The coupling between nearby superconductors sepa- 
rated by a ferromagnet is a property that follows from 
the proximity effects, which in the context of F/S multi- 
layers consist of the existence of superconducting correla- 
tions in the ferromagnet and magnetic correlations in the 
superconductor, arising from their mutual influence. The 
resulting superconducting phase coherence is quantified 
by the pair amplitude F(r) — (tpi(r)ij}^(r)), where the 
ij)„ are the usual annihilation operators. It is now well 
established that the leakage of superconductivity is due 
to the process of Andreev reflection^, whereby a quasi- 



particle incident on a F/S interface is retroreflected as 
a quasihole of opposite spin. It is in turn the coherent 
superposition of these states, spin split by the exchange 
field in the ferromagnet, that ultimately leads to damped 
oscillations of F(r) in the magnet, with a characteristic 
length £p typically much smaller than the superconduct- 
ing coherence length £o- These oscillations arc akin to 
high field oscillatory phenomena described a long time 
agoAi In the absence of currents and magnetic fields, the 
modulation of the order parameter determines whether 
two neighboring superconductor layers share a stable tt 
or phase difference. For a multilayer F/S heterostruc- 
ture with ferromagnet layers of order tt^f in width, it is 
intuitively evident, taking into account the continuity of 
F(r) across the F/S interface and the particular oscilla- 
tory nature of the pair amplitude in the ferromagnet, that 
a configuration will result in which it is energetically fa- 
vorable to have a phase difference of A</> = tt, rather than 
zero, between successive superconducting layers. This in- 
deed turns out to be the case. 

Although recently there has been a surge of inter- 
est in the study of F/S multilay er str uctu res, (se e e.g 
the theoretical work of Refs. IS U9llOUlllll2lll3Ul4lll5lllt: 
[T7lfl^fT9ll20ll2]]|23l2M2^l25l and the experimental work 
discussed below) work on superconductor-ferromagnet- 
superconductor (SFS) Joscphson junctions started long 
agoS The Josephson current was calculated for a short 
weak link in the clean limit, and found to exhibit oscilla- 
tions as a function of the ferromagnet exchange fieldiSi It 
was later demonstrated that for an SFS sandwich obey- 
ing the dirty limit conditions, the critical current os- 
cillates as a function of the thickness of the magnet, 
and of the exchange fields A more detailed analysis 
of dirty tt junctions near the critical temperature al- 
lowed for differing transparencies of the ferromagnet- 
superconductor interfaces^ Many interesting phenom- 
ena have been proposed or discussed. Calculations were 
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more recently performed for an SFS junction with ar- 
bitrary impurity concentration. For non-homogeneous 
magnetizationji£i2& the superconductor may exhibit a 
nonzero triplet component extending well into the mag- 
net. For quasi two-dimensional, tight-binding, F/S 
atomic-scale multilayers, the ground state was shown 
in some cases to be the tt stated, and the density 
of states (DOS) exhibited prominent features that de- 
pend critically on the exchange field and transfer inte- 
gral parameters^ For multilayer structures consisting of 
two ferromagnets and an insulator sandwiched between 
two superconductors, an enhancement of the Josephson 
current was predicted22i2i for antiparallel alignment of 
the magnetization in the ferromagnet layers. Spin-orbit 
scatterin g 1 ^' 22 and changing the relative orientation an- 
gle of the in-plane magnetizations^ 3 , were shown to signif- 
icantly modify the behavior of the dc Josephson current. 

The rapidly evolving theoretical views compounded 
with technological advancements which permit the 
fabrication of well-characterized heterostructures, has 
prompted a considerable number of experimental inves- 
tigations of 7r coupling on several fronts. A study of the 
superconducting transition temperature for F/S multi- 
layers revealed oscillatory behavior as a function of fer- 
romagnet thickness*^ For tt junctions involving relatively 
weak ferromagnets, variations in temperature can induce 
a crossover from Atfi — to the tt state, and this was 
observed^ as oscillations of the critical current versus 
temperature. The transition to the tt state is also re- 
flected in critical current measurements for Josephson 
junctions in which the ferromagnet layer separating the 
two superconductors was systematically varied. 30 The 
superconducting phase was measured directly^ using 
SQUID's made of tt junctions, demonstrating a half quan- 
tum flux shift in the diffraction pattern. Direct evidence 
of the oscillatory behavior of the superconducting corre- 
lations in the ferromagnet was found through tunneling 
spectroscopy measurements which yielded inversions in 
DOS for a thin ferromagnetic film, in contrast with the 
behavior in a superconductor^ 

A common feature that pervades most of the theo- 
retical work mentioned above is the use of quasiclas- 
sical formalisms, often compounded by the neglect of 
self-consistency for the space dependent pair potential, 
A(r). These approximations do have the advantage of 
providing an accessible and efficient method to approxi- 
mately calculate properties of inhomogeneous supercon- 
ducting systems, while avoiding the cumbersome numer- 
ical issues that arise when attempting to solve the cor- 
responding, much more complicated, self-consistent mi- 
croscopic equations. The general underlying drawback 
of such approximations however, is the elimination from 
consideration of phenomena at the atomic length scale 
given by the Fermi wavelength, Xp, as can be seen in 
the derivation of the Eilenberger equations. 33 Further ap- 
proximations follow when the assumption is made that 
the mean free path is much shorter than £07 m which case 
the Eilenberger equations reduce to the widely used Us- 



adel equations! 3 * The elimination of the relatively small 
length scales poses problems for quasiclassical methods 
(even when self-consistent) when interfacial scattering 
is involved^ or when the geometry or potentials have 
sharp variations on the atomic scale. These issues, of 
increasing experimental importance given the ever im- 
proving quality of the experimental samples, often re- 
quire nontrivial effective boundary conditions that must 
supplement the basic equations. The problem worsens 
when dealing with multilayer structures, where the suc- 
cessive reflections and transmission of quasiparticles cre- 
ates closed trajectories 3 ^ that may render the quasiclas- 
sical approximation scheme inapplicable. 

In this paper we investigate the proximity effect and as- 
sociated electronic properties of clean three-dimensional 
F/S multilayer structures comprised of alternating super- 
conductor and ferromagnet layers. Our emphasis is on 
the study of the existence of pair potential behavior of 
the tt type. We implement a complete self-consistent mi- 
croscopic theory that treats all the characteristic length 
scales on an equal footing, and thus can accommodate 
all quantum interference effects that are likely to be per- 
tinent. The problem will be solved from a wave function 
approach using the Bogoliubov-de Gennes (BdG) equa- 
tions. To do so, we extend an earlier method 3 ^ used 
for a single F/S structure, to allow for a more compli- 
cated geometry, consisting of an arbitrary number of lay- 
ers. Self-consistency is rigorously included, as it has been 
demonstrated 3 ^ 3 ^ that this is essential in the study of the 
proximity effect at F/S interfaces. 

We present in Sec|H]the geometry and the numerical 
approach we take to solve the microscopic BdG equations 
and obtain the self consistent energy spectra (eigenval- 
ues and eigenfunctions). We also explain in some detail 
how the local DOS and the ground state energy are cal- 
culated. In Sec lIIII we first present our results for SFS 
structures and show examples of the relevant quantities: 
these include first of all the pair amplitude, which is used 
to illustrate the cases in which either the or tt state is 
energetically favored. The thickness of the F layer or 
layers turns out to be the decisive parameter, with the 
zero and tt states periodically alternating in stability as 
this quantity varies. The experimentally accessible DOS 
averaged over a superconducting layer is next discussed: 
results for both the sum and the difference of the up and 
down spin terms are presented and their correlation with 
the zero or tt states demonstrated. Results for the local 
magnetic moment, which we show how to calculate from 
the spin-dependent local DOS, are also given: in the su- 
perconductor this quantity measures the penetration of 
magnetic correlations. We analyze also, in a similar fash- 
ion, a more complicated five layer structure, discuss the 
similarities and differences between the two geometries, 
and the generalization of our results to more complicated 
structures. Finally in Sec II VI we briefly summarize our re- 
sults and discuss potential experimental implications and 
future work. 
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FIG. 1: Schematic of the model geometry used in this paper. 
The total thickness in the z direction is d, and the thicknesses 
of the S and F layers are ds and dF as indicated. There is a 
total number Nl of layers Nl = 3 and Nl — 5 in this work, 
with the outer ones being superconducting. The line breaks 
in the middle region denote repetition. 



II. METHOD 

In this paper we consider a semi-infinite multilayer 
structure of total length d in the z-direction, consisting 
of an odd number Nl of alternate superconductor (S) 
and ferromagnetic (F) layers, each of width ds and d F 
respectively (see Fig. QJ. The sandwich configuration is 
such that the complete structure begins and ends with 
a superconductor layer. The free surfaces at z = and 
z = d are specularly reflecting. The basic methodology 
we use is an extension of that which has been previously 
discussedi^SiSi Upon taking into account the translational 
invariance in the x — y plane, one can immediately write 
down the BdG 38 equations for the spin-up and spin-down 
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where the free-particle Hamiltonian is defined as, 

E F (z). (2) 



2m dz 2 



Here e± is the transverse kinetic energy, the e„ are the 
quasiparticle energy eigenvalues, and A(z) is the pair po- 



tential, described below. The magnetic exchange energy 
ho(z) is equal to a constant ho in the ferromagnet lay- 
ers, and zero elsewhere. A potential U(z) describing in- 
terface scattering can easily be added to Eqn. [3 We 
define the quantity Ep{z) to equal Efm m the mag- 
netic layers, so that in these regions, Ep\ — Efm + ho, 
and Efi = Efm — ho- Likewise, in the superconduct- 
ing layers, Ep(z) — Eps- The dimensionless parameter 
/ = ho/ Efm characterizes the strength of the magnet. 
At / = 1, one therefore reaches the half metallic limit. 
From the symmetry of the problem, the solutions for the 
other set of wavefunctions (u^v^) are easily obtained 
from those of Eqns. (|TJ by allowing for both positive and 
negative energies. The BdG equations are completed by 
the self consistency condition for the pair potential, 



A(z) 



gi z ) 

2 



[ul(z)vi(z) + ui(z)vl(z)] tanh(e„/2T), 

(3) 



where T is the temperature, g(z) is the effective cou- 
pling describing the electron-electron interaction, which 
be take to be a constant g within the superconductor lay- 
ers and zero within the ferromagnet layers, and ojd is the 
Debye energy. We have not included spin-orbit coupling, 
and assumed that all of the F layers are magnetically 
aligned and hence2^ considered singlet pairing only, in 
the s-wave. 

We solved Eq. by expanding the quasiparticle am- 
plitudes in terms of a finite subset of a set of orthonor- 
mal basis vectors, u\(z) = ^2 q u} iq 4) q (z), and v*(z) = 
12 q v nq ( l > q( z )- W e use the complete set of eigenfunctions 
(j)q(z) = (z\q) = y/2/dsm(k q z), where k q = q/nd, and q 
is a positive integer. The finite range of the pairing in- 
teraction Lop, permits the number N of such basis vectors 
to be cut off in the usual way^ Once this is done, we 
arrive at the following 27V x 2N matrix eigensystem, 



H+ D 
D H~ 



(4) 



where ^ 



' ,- 1 • • • • i U nN' 

ements H*, connecting <j) q to 4> g > 
the real-space quantities in Eq.JTJ, 



, v^ N ). The matrix el- 
are constructed from 



J_d 2 

2m dz 



ho{z) 



k 2 
2m 



2 , ej_ - E F (z) 

d r d 



S qq > - dz(j>q(z)E F i(z)(f> q i(z) - / dz4> q (z)E F {z)4> q <(z). 



(5a) 



The expression for H , is calculated similarly. The off-diagonal matrix elements D qq i are given as, 



D gg , = (q\A(z)\q')= / dz<p q (z)A{z)<t> q ,(z). 
Jo 



(5b) 
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After performing the integrations, Eq.JSaJ can be expressed as 

1 \ E F y r sin[n(fc 9 - k q >)(d F + d s )/2] sm[(k q - k q >)(n(d F + d s )/2 - d F )} 



[kg k q r) 



[k q k q i ) 



3in[(fc g + k q ,){n{d F + d s )/2 - d F )} sin[n(fc g + k q >)(d F + d s )/2] 



E F s 



[kq + kqi) (kq + kqi) 

sin[(fc g — k q ')(n(d F + d s )/2 - d F )] sin[n(A; g - k q i)(d F + ds)/2] 



(k q k q ') 



(kg k q i) 



sin[n(fc g 4- k q i){d F + d s )/2] sin[(fc 9 4- k q >)(n(d F + d s )/2 - d F )] 



(k q + kqi) (k q + k q i) 

where the sum is over even integers only. The diagonal matrix elements are somewhat simpler, and are written, 

E F ^ 



(6a) 
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N L -1 



NLd F — — sm[nk q (d F + ds) — d F ] — sm[nk q (d F + dg)] 



E 



FS 



2,1 



N L -1 



NLds — — E sm[nk q (d F + d s )] — sm[nk q (d F + ds) - d F ] 



The self-consistency condition, Eq.Q, is now transformed into, 

tt\(z) 



A(z) 



k F d 



p,p' q 



d£± u "P V np' + u n P v lp> sin(fc p z)sin(Vz)tanh(e n /2T), 



(6b) 



(7) 



where \(z) = g(z)N(0), and N(0) is the DOS for both 
spins of the superconductor in the normal state. The 
quantum numbers n encompass the continuous trans- 
verse energy s± , and the quantized longitudinal momen- 
tum index q. 

The primary quantity of interest is the local density of 
one particle excitations in the system, N(z,e). Current 
experimental tools such as the scanning tunneling mi- 
croscope (STM) have atomic scale resolution, and make 
this quantity experimentally accessible. Since it is as- 
sumed that well defined quasiparticles exist, the tun- 
neling current is simply expressed as a convolution of 
the one-particle spectral function of the STM tip with 
the spectral function for the ferromagnet-superconductor 
system^ The resultant tunneling conductance, which is 
proportional to the DOS, is then given as a sum of the 
individual contributions to the DOS from each spin chan- 
nel. We have: 

N(z,s) = N 1 (z,e) + N l (z,e), (8) 
where the local DOS for each spin state is given by 

N r (z,e) =-^{[^)] 2 /'(e - e») + [vl(z)] 2 f'(e + e n )}, 

n 

(9a) 

e) = - J2{[ui(z)ff'(e - e n ) + [vi(z)] 2 f'(e + e„)}. 

n 

(9b) 



Here thermal broadening is accounted for in the term 
involving the derivative of the Fermi function /, /'(e) = 
df/de. 

We shall see below that we will also need to compare 
different self-consistent states. In general this is done in 
terms of the free energy. However, we will consider here 
only the low temperature limit. For T — > 0, the entropy 
term can be neglected, as it vanishes proportionally to 
T 2 . In this case all we need is the ground state energy 
Eq . In evaluating this quantity some care must be taken 
in properly including all energy shifts, even in the bulk 
In the inhomogeneous case the result^ can be 



eN(z,e)de+-(\A(z)\ 2 ) (10) 
9 
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cas' 
written as: 



E Q = dz 
Jo 



where the angular brackets in (| A(z)| 2 ) denote the spatial 
average, and N(z, e) is given in Eqn. [S] One can rewrite 
Eq in a somewhat more standard way: 



E Q 



;(|A(z)| 2 ), (11) 



which in principle gives E in terms of the calculated 
excitation spectra. 
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III. RESULTS 

In this section we present and discuss the results that 
we have obtained through our numerical solution of the 
matrix eigensystem Eq. Q and the self-consistency con- 
dition Eq. Q) . We will study the two cases of Nl = 3 and 
Nl = 5, that is, SFS and SFSFS structures separately. 
We consider only "regular" structures in which all S lay- 
ers have the same thickness ds, which we will take to be a 
fixed value larger than £o 5 while the F layers, when there 
is more than one, have all the same thickness, dp, which 
we will vary. With the assumption that no current flows 
across the sample, the quantity A(z) can be taken to be 
real, but it can in principle switch sign (zero or tt state) 
in going from one S layer to the next. 

In our calculations we have studied two different values 
of the parameter /, / = 0.5 and I — 1. We have set the 
superconducting correlation length £o to So = fcs£o = 50, 
where k$ is the Fermi wavevector of the superconduc- 
tor, and taken lu = ujo/Eps = 0.1 for the dimensionless 
Debye energy cutoff. It follows from previous studies^ 
that the first of these parameters simply sets the over- 
all length scale in the superconductor and is of little 
relevance whenever dg exceeds £o ; as will be the case 
here, while the second is unimportant at low tempera- 
tures (the limit that we will consider), as it simply sets 
the scale for T c . We have also assumed that there is no 
oxide barrier between the layers and that the "mismatch 
parameter" A = (Efm I 'Eps) is unity A nonzero bar- 
rier height would in general diminish the amplitude of 
all the phenomena discussed here, without qualitatively 
altering the results. The possible influence of varying A 
is more complicated: this parameter^ determines, to- 
gether with /, the basic spatial periodicity of the prob- 
lem, £p k (fef — /c|) _1 , (where fef and k± are respectively 
the Fermi wavevectors of the parabolic up and down spin 
bands in the ferromagnet) which we shall see is very im- 
portant here. Furthermore, the amplitude of the oscilla- 
tory behavior found in simpler SF structures decreases^ 
with A. 

As explained in previous work (see Refs. the 
self consistent solution to these equations is obtained it- 
eratively: one makes a suitable initial guess for A(z), di- 
agonalizes the system Eqn. (@J for that guess, and com- 
putes an iterated A(z) from 0). The process is then 
repeated until convergence is obtained. The technical- 
ities for the self consistent solution of these equations 
were extensively discussed in previous work^Lii. The di- 
agonalization in terms of the orthonormal basis chosen 
must be performed for each value of e± in the appropri- 
ate range. We took here N± = 5000 different values of 
e± , except as indicated below, and the number of basis 
functions required for convergence was up to N = 1000. 
The self consistent solution process is terminated when 
the relative error between consecutive iterated values of 
A(z) nowhere exceeds 10~ 4 . We have found that the 
number of iterations needed to achieve self consistency 
can be quite large: in most cases, it exceeds fifty. 



Because our objective here is to discuss the possible ir 
states, in starting the iteration process we make two dif- 
ferent initial guesses: one is of the ordinary "zero" state 
form, where the initial guess has the same sign (conven- 
tionally positive) in all the superconducting layers, and 
one of the tt form, where it alternates sign from one S 
layer to the next. We have found that in some cases, for 
example for SFS structures with small ds (i.e., < £o)j and 
dp < 7t£f, the self-consistent A(z) typically converges to 
either a or 7r-state regardless of the initial guess, de- 
pending on (if. A similar trend holds for the small ds 
five layer SFSFS system but over a broader dp range. 
However, for the regular structures that we will focus on 
here, with ds 3> £o> two different self consistent solutions 
are always obtained, one of the zero and one of the tt 
type, according to the type of initial guess. We interpret 
this as showing that two local minima of the free energy 
exist. We then have to determine the stable minimum 
by calculating the free energy (or rather, at low temper- 
ature, the ground state energy) of both self consistent 
states, as discussed below, and comparing them. 



A. SFS 

We consider first the case of an SFS sandwich. Prelim- 
inary investigations showed that the situation of interest 
occurs when the magnetic layer is not too thick. This 
is as expected, since the overall length over which the 
superconducting correlations penetrate (in an oscillatory 
way) into the magnet is characterized by (fcj — k^ 1 , 
which at the relatively large values of I considered here 
is fairly small. Thus, we have taken in the studies pre- 
sented here a thickness ksds = 300 for the superconduct- 
ing layers, and the parameter ksdp meanwhile, is varied 
in the range between one and twenty. The choice of ksds 
determines, through standard BCS theory relations, the 
value of the ratio of the superconductor Fermi energy to 
the bulk order parameter. 

As explained above, results were obtained by iteration 
from two initial configurations of A(z), with the initial 
guesses corresponding to opposite signs for the pair po- 
tential in each of the two S layers. For the range of 
parameters considered here, both initial guesses led in 
all cases to self consistent configurations, which were ei- 
ther of the zero or of the tt types, according to the initial 
guess. This is described in Fig. [21 where we show ex- 
amples of the two self consistent solutions for the pair 
amplitude F(Z), as a function of the dimensionless dis- 
tance Z = ksz. It is not surprising that both types of 
solutions are found: it is after all obvious that in the limit 
where dp is sufficiently large, both solutions must exist 
and be degenerate. The pair amplitude F(z) = gA.(z) 
does not vanish identically in the magnetic region, but 
it exhibits the well-known oscillations. In the supercon- 
ductor, it rises in absolute value towards the bulk result, 
away from the S/F interfaces. Results are shown for two 
values of I and two values of ksdp. We can see that in 
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FIG. 2: (Color online). Results for the pair amplitude F(Z), 
normalized to the bulk superconductor value, as a function of 
Z = ksz, in an SFS structure. The dimensionless thickness of 
the S portions is ksds = 300, while the corresponding values 
of the dimensionless thickness of the intervening F layer are 
(from top to bottom) ksd,F = 10, 16, 19. The blue (solid) lines 
represent self-consistent solutions of the zero type, and the red 
lines alternative self consistent solutions of the tt type. The 
value of I is 0.5 and the superconducting correlation length 
is £ = 50. 



certain cases, depending on the thickness of the F layer, 
F(z) in the superconductor (and hence A(z)) is larger, 
in absolute value, for the 'zero' than for the tt state (see 
top panel), while in some other cases (middle panel) the 
opposite occurs, and for yet some other dp values (see 
the bottom panel) there is no observable difference. In- 
tuitively, this happens because of the different way, de- 
pending on ksd,F, in which the pair amplitude in the two 
superconductor regions must adjust itself to the oscilla- 
tions in the magnet. That the oscillatory behavior of the 
pair amplitude in the F layer is clearly different for the 
zero and tt solutions can be seen by careful examination 
of the portion of the plots which lies in the F region. 

To find out the most stable configuration, one must 
compute the difference in the ground state energies, or 
equivalently the condensation energies, of the zero and 
tt states. This can be done in principle by using Eqn. 
(|ll|l . In practice this is computationally very difficult: 
the value of Eq for each state must be computed sepa- 



rately from its own spectrum (which can consist of up to 
10 6 eigenstates) , and the results subtracted. Since each 
Eq includes the normal state energy, which is many or- 
ders of magnitude larger than the condensation energy 
sought, this requires extreme numerical accuracy. The 
problem is exacerbated because the "logarithmic" last 
term in the right side of Eqn. I|ll|) is in itself much 
larger than the condensation energy (and the latter is 
itself considerably larger, as we shall see below, than the 
condensation energy difference between the two states), 
and must be exactly canceled by a portion of the first 
term. This is a well-known problem even in the bulk 
case, where great care has to be takenffi to make the del- 
icate cancellation analytically explicit. We have found it 
technically impractical to numerically compute Eq from 
Eqn. I|ll|) for all cases considered^ with the required 
precision. However, by using increased values of N± and 
N in a few selected cases we have been able to verify that 
the ground state condensation energy (that is, after sub- 
tracting the normal ground state energy E 0n calculated 
for the same geometry and parameter values except for 
setting g — 0) for either the zero or tt states is, for the 
cases considered here where dp is small and d$ >> £o, 
approximately given by: 

E Q -E Qn ^-aN(0)(\A\ 2 ). (12) 

This result is, a posteriori not surprising at all in the limit 
of large c?s and small dp, as it is quite similar to what is 
found^ analytically for the bulk: in that case a is exactly 
0.5 and the spatial average is of course replaced by the 
uniform bulk value. In our case we find the coefficient 
a < 0.5 within our numerical uncertainty. The right side 
of Eqn. H12JI is of course very easy to compute. Thus, we 
have adopted a procedure based on Eqn. i|12|l to compare 
condensation energies for the two competing states. 

The results are shown in Figure |31 The quantity plot- 
ted there is the difference between the values of (|A| 2 ) for 
the zero and tt states normalized to iV(0)A§, where Ao 
is the bulk gap. This normalization corresponds to twice 
the bulk value limit of the condensation energy. This 
is then a dimensionless measure of the condensation (or 
equivalently, ground state) energy difference between the 
self consistent zero and tt configurations, (see Eqn. dJ). 
This normalized energy difference is plotted as a function 
of the dimensionless thickness ksdp of the intermediate 
F layer, which is sandwiched between thick (ksds — 300) 
S layers. We see that the difference in energies is, as one 
would expect, only a small fraction (about one tenth at 
the most) of the bulk condensation energy. We also see 
that it is an oscillatory function of ksdp. Comparison of 
the top and bottom panels (which correspond to I = Q.5 
and 1 = 1 respectively) shows that the rough periodicity 
of these results is approximately given by (fcj — fcj) -1 , and 
it is in fact very similar— ^ quantitatively to the oscillatory 
behavior of the pair amplitude F (z) in a thick magnetic 
layer. At small ksdp, the zero state is obviously very 
favored, as one would expect, while in the limit of large 
ksdp the energy difference is of course zero, reflecting 
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FIG. 3: (Color online). The difference in condensation ener- 
gies, AEo, between the zero and tt states for an SFS sandwich 
in the low temperature limit, normalized to N(0)Aq, calcu- 
lated as explained in the text. The results are plotted as a 
function of the dimensionless thickness ksd,F of the ferromag- 
netic layer, for two values of /. At small cLf the zero state 
is favored. The periodicity of the results is determined by 
(fcj — fej.) -1 , as expected. 



the degeneracy of the two states. The influence of the 
parameter / is quite dramatic: in the half metallic case 
(lower panel) the first peak favoring the 7r state is more 
prominent and the zero state is generally speaking less 
favorable than for the intermediate value of / shown in 
the top panel. 

We turn now to the density of states (DOS) for this ge- 
ometry. Typical results are exhibited in Fig. 01 where we 
show the DOS, integrated over the superconducting re- 
gion of thickness ksds = 300, as a function of the energy 
normalized to Ao- Results are shown for two values of / 
(top and bottom panels) and, for each value of /, at two 
values of ksdp, one corresponding to the case where the 
equilibrium state is of the zero type, and the other cor- 
responding to the opposite situation. One can see that 
for / = 0.5 there are states in the gap, and that these 
states are more prominent in the 7r case where there is 
a zero energy small peak. At / = 1 the DOS results are 
also different: although for both of the cases shown there 
is a gap in the spectrum, the location of the peaks near 




FIG. 4: (Color online). Density of states (DOS) results for 
SFS structures. The quantity plotted is the local DOS inte- 
grated over one S layer, normalized to iV(0). The energy is 
normalized to the bulk gap Ao. The top panel shows results 
at ksd,F = 5, where the stable state (see Fig. |3J| is of the tt 
type (red curve, labeled as tt) and at ksd,F = 10, where the 
zero state is more stable (blue solid curve, labeled "0"). In 
the bottom panel, 1 = 1 and consistent with the doubling of 
/, the thicknesses displayed are halved to ksd-F = 2.5 (it case) 
and ksd-F = 5 (zero case). See text for discussion. 



the gap edge is not the same for the zero and tt states, 
with the first peaks being more prominent and at higher 
energies in the latter case. Thus, there are genuine differ- 
ences between the DOS of zero and tt states, which may 
be experimentally observable. 

It is also of interest to show the difference between the 
local DOS for up and down states, as defined by 



5N(z,s) = N 1 (z,e)-N l (z 1 e). 



(13) 



This is done in Fig. [5] where results are shown for the 
two cases corresponding to those also displayed in the top 
panel of Fig.0J The differential DOS shown is integrated 
over the thickness of one S layer, and normalized to the 
total normal bulk DOS value. Because of the finite value 
of /, the results are not symmetric around zero energy. 
One can see that the energy structure at the gap edge is 
appreciably more prominent for the thickness value which 
corresponds to an equilibrium tt state, while for the zero 
state the structure is broader and more diffused. 

An alternative way of illustrating the magnetic polar- 
ization effects, which has also the advantage of providing 
local information, is through the use of the local mag- 
netic moment m(z). This quantity is easily obtained by 
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FIG. 5: (Color online). Differential density of states (DOS) 
between up and down spin states for SFS structures. The 
quantity plotted is that defined in Eq. ilty . integrated over 
one S layer, and normalized to N(0). Results are shown for 
/ = 0.5 at ksd,F = 5, where the stable state (see Fig. [3J is 
of the 7r type (red curve) and at fes^F = 10, where the zero 
state is more stable (blue solid curve). 

integration of the local DOS results. One has: 

m(z) = fi B J deSN{z,s)f(e), (14) 

where /is is the Bohr magneton and the integral extends 
over the occupied states in the band. This can be cast in 
a more convenient form as: 

m(z)= Mb^tO)) - (ni(z))], (15) 

where (n a (z)) is the average number density for each 
spin subband, and is written in terms of the quasiparticle 
amplitudes as, 

(M*)> = E{K(^)] 2 /(e»)+KW] 2 [l-/(e„)]}, a =T 

n 

It is more instructive to plot m(z) normalized to the cor- 
responding integral of N^(z, e)+N±(z, e). We denote this 
normalized quantity by M(z) and we plot it, in units of 
the Bohr magneton, in Fig. EI The two panels there cor- 
respond to values of I and ksd,F as in the corresponding 
top and bottom panels of Fig. 0] We see in this figure 
that for relatively small ksdp, the quantity plotted rises 
up sharply from the F/S interface and then has a slow 
modulation as it approaches its bulk value in the F layer. 
The magnetization does not vanish identically inside the 
superconductor: its behavior there consists of strongly 
damped oscillations, with an overall characteristic spa- 
tial decay on the order of a few Fermi wavelengths. The 
effect does not seem to depend strongly on whether one 
is dealing with zero or it states. 

The self consistent results displayed can also be inter- 
preted as representing an effective, local value of I{z), 
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FIG. 6: (Color online). Normalized local magnetic moment 
as defined in the text and Eqn. 1151 Results in the top and 
bottom panels correspond to the same values of I and thick- 
ness as in the corresponding panels of Fig. [I] Thus the top 
panel is for I = 0.5 and ksd-F = 5, 10, while the bottom panel 
is for I = 1 and ksdp = 2.5, 5. 

through the relation M(z) = /j, B [(l + /(z)) 3 / 2 - (1 - 
I{z)fl 2 ]/[{\ + I{z)f' 2 + (1 - I(z)f/ 2 }. The quantity 
I(z) is then the magnetic counterpart of the self consis- 
tent F(z), measuring directly the magnetic part of the 
proximity effect, that is, the leakage of magnetic correla- 
tions into the superconductor. 

B. SFSFS 

: 1 • 

In this subsection we consider the case of more com- 
plicated, five layer structures. These are realizable 
experimentally^ and therefore of considerable interest. 
As in the three layer case, we will study the situation 
where the three superconducting layers are relatively 
thick, taking again ksds = 300 and the F layers are 
thin enough so that F/S proximity effects cannot be ne- 
glected. 

We begin by considering (see Fig.EJ) the pair amplitude 
F(z). This figure is in every way analogous to Fig- El ex_ 
cept for the insets, where we display in more detail the 
behavior of F(z) in one of the ferromagnetic layers. Re- 
sults for solutions of both the zero and the 7r type are 
shown. Both are obtained self-consistently, the first by 
starting from an initial guess in which the sign of the or- 
der parameter in the three S layers is always the same, 
and the second by starting with a guess in which the or- 
der parameter in the middle S layer is the opposite to 
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FIG. 7: (Color online). Results for the pair amplitude F(Z), 
normalized to the bulk superconductor value, as a function of 
Z = ksz, in an SFSFS structure, for I = 0.5. The dimen- 
sionless thickness of the S portions is ksds = 300, and the 
corresponding values of the dimensionless thickness of the in- 
tervening F layer are (from top to bottom) ksdF = 10, 16, 19. 
The blue (solid) lines represent self-consistent solutions of the 
zero type, and the red lines of the 7r type. The insets are a 
magnification of one of the F regions. The vertical axis in the 
insets varies between ±.35 in dimensionless units. 



that in the other two layers. Self consistent solutions are 
always reached by iteration, for large ds and cLf in the 
ranges shown, in either case. We observe the expected 
depletion of F (z) near the F/S interfaces, and the sub- 
sequent approach towards it bulk value over the length 
scale £o> with the maximum A(z) in the central S layer 
slightly reduced from the bulk Ao. We also see in the 
main panels that depending on the value of kgdp, the 
absolute value of F(z) in the superconductors varies pe- 
riodically between being larger in the zero state to be- 
ing larger in the it state, as was the case for three layer 
structures. The insets illustrate more clearly how the ex- 
istence of the two states relates to the oscillations of F{z) 
in the ferromagnetic region, which are very different in 
each case. 

As in the three layer case, therefore, we find that there 
are two local minima of the free energy, corresponding 
to the zero and 7r alternatives. Again, the absolute min- 
imum, at low temperature, must be found by comparing 



FIG. 8: (Color online). Difference in condensation energies, 
AEo, between the zero and n states for a five layer SFSFS 
system, calculated as explained in the text, and normalized 
to jV(O)Aj), as in Fig.|H] This quantity is plotted as a function 
of the thickness ksdF of each ferromagnetic layer. Results for 
two values of I are shown. 



the two condensation energies. This we do in the same 
way as for the three layer case (see Eqn. ^] and associ- 
ated discussion). The results are shown in Fig. [8] which 
should be compared with Fig. |3J The two figures are 
remarkably similar. In both cases the behavior is os- 
cillatory, with the same approximate spatial periodicity 
related to that of the pair amplitude oscillations. Again, 
the obvious results that the zero state is favored at small 
k$dp and that the two states are degenerate for large 
ksdF are recovered. The three and five layer plots are 
not identical, however: in the latter case we find that 
the overall scale of the phenomenon is nearly a factor 
of two higher, as one can see by comparing the vertical 
axes. The first peak favoring the ir state is higher and 
sharper for five layers. Although the effect of increasing 
the layer number is not as dramatic as that of increasing 
I, one can nevertheless assert from the trend that the 
oscillatory behavior with dp would not only persist but 
would be even more prominent if the number of layers 
were further increased, as in superlattices. 

A few selected DOS results for this SFSFS geometry 
are shown in Fig. El which should be viewed in compari- 
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FIG. 9: (Color online). Density of states (DOS) results for 
SFSFS structures. The local DOS integrated over one of the 
external S layers, normalized to N(0), is plotted vs. the en- 
ergy normalized to the bulk gap Ao. The top panel shows 
results at ksdF = 5, where the stable state (see Fig. |SJ is 
of the 7r type (red curve, labeled as 'V") and at ksdF — 10, 
where the zero state is more stable (blue solid curve, labeled 
"0"). In the bottom panel, 1 = 1 and the thicknesses are 
ksdF = 2.5 (n case) and ksdF = 5 (zero case). 

son with the analogous Fig.0]for the SFS structure. The 
quantity plotted is averaged over one of the two outside 
S layers, and all parameters are chosen to be the same 
as in Fig. 0] The similarity between the two figures is at 
first sight very remarkable, although a second look shows 
that the structure of the subgap peaks is far from being 
the same, particularly for the zero state case, where ad- 
ditional shoulders appear at / = 1. One concludes again 
that many features, including the zero energy peak in the 
stable 7r state of the top panel, are robust with respect to 
increasing the number of layers, and very likely to persist, 
and even be more obvious, in larger regular structures. 

The differential DOS between up and down states for 
this geometry exhibits a behavior sufficiently similar to 
that displayed in Fig. for the SFS case that there is no 
need to display it in a separate figure here. On the other 
hand, it is worthwhile to illustrate an example of the 
normalized local magnetic moment M{z). This is done 
on Fig. where this quantity, as defined in Eqn. ITU 
is plotted with the same normalization and parameter 
values as in the top panel of Fig. The behavior for 
the two geometries is certainly similar, but one again 
sees that the magnetic penetration effects become more 
prominent as the number of layers increases from three 
to five. This is another indication that such effects are 
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FIG. 10: (Color online). Normalized local magnetic moment 
for an SFSFS structure, as defined in the text and Eqn. 1141 
Results are for / = 0.5 and ksdF = 5, 10. 



very likely to be easier to observe in structures involving 
a larger number of layers. 



IV. CONCLUSIONS 

We have rigorously investigated the proximity effects 
that occur in clean multilayered F/S structures of the 
SFS and SFSFS type. We used a microscopic wave- 
function approach that does not coarse grain over length 
scales of order A f , and thus accounts for atomic-scale ef- 
fects. The space dependence of the pair amplitude F(z) 
was obtained self-consistently by using an efficient nu- 
merical algorithm. From the calculated eigenstates, we 
were then able to obtain the experimentally relevant local 
magnetic moment, and the local density of states. 

We have demonstrated that for all the cases considered, 
where the thickness ds of the superconducting layers is 
much greater than £o and that of the ferromagnetic re- 
gions is relatively small, two local minima of the ground 
state energy exist, thus yielding self-consistent states of 
the and n types. Through a careful analysis of the pair 
amplitude and excitation spectrum, we have calculated 
which of these two states is the actual ground state, with 
the lowest energy. The results show that the difference in 
condensation energies between the and tt states exhibits 
damped oscillations as a function of ferromagnet width, 
with the characteristic exchange-field dependent spatial 
period being given approximately as 27r(fc-f — the 
same quantity which characterizes the oscillations of F{z) 
in bilayers. The local DOS exhibits strikingly different 
behavior for two exchange fields that differed by a fac- 
tor of two. For I = 0.5, the subgap DOS shows a gapless 
structure, with features that depend strongly on whether 
the ferromagnet width corresponds to the or 7r state. 
The half-metallic case (I = 1.0) is on the other hand 
gapless in the range of dp considered, and the modified 
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excitation spectrum reveals itself through the differing 
peaks in the DOS. To illustrate the leakage of magnetism 
into the superconductor, the differential DOS between 
the spin up and spin down states was presented for a 
SFS junction. The most prominent spin-splitting was 
seen for the 7r-junction at energies e/Ao ~ 1. We believe 
that this represents an experimentally important signa- 
ture for the 7r state. We have also calculated the local 
magnetic moment for both the three and five layer cases, 
to give further insight into magnetic polarization effects. 
Although we found the results to be relatively insensitive 
to a or 7r state configuration, we were able to extract 
an effective local value of I(z) in both the F and S layers. 

The calculations and method used in this paper, al- 
though sufficiently general to include in the future more 
complicated effects (e.g. finite temperature, other pairing 
states, spin- flip scattering, and impurities), were taken 
within the ballistic limit. This limit is appropriate for 
ferromagnet layers whose width is less than the mean free 



path, and this is consistent with our calculations, where 
we have taken kgdp < 20. The inclusion of interfacial 
scattering would likely have the effect of diminishing the 
proximity effect, without qualitatively altering the char- 
acteristic results. For bulk impurity scattering, the ef- 
fective £f would involve not just / but also the diffusion 
length. It is the goal of future work to address these top- 
ics, and also others, including heterostructures comprised 
of a single superconductor sandwiched between two ferro- 
magnets with arbitrary relative magnetization, F/S mul- 
tilayers with a greater number of layers, and smaller su- 
perconductor widths, where geometrical and atomic-scale 
effects are likely to be more prevalent. 
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